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Abstract 

A free-parameter linear-combination-of-atomic-orbitals approach is presented 

for analyzing the stopping power of slow ions moving in a metal. The method 

is applied to the case of He moving in alkali metals. Mean stopping powers for 

He present a good agreement with local density approximation calculations. 

Our results show important variations in the stopping power of channeled 

atoms with respect to their mean values. 
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I. INTRODUCTION 



The interaction of ions with condensed matter has drawn the attention of many re- 
searchers from the beginning of this century |l[ . A great deal of work in this field has dealt 
with the energy loss of swift ions in solids. In this regard the work of Bethe , Fermi || , 
Williams H, and Lindhard || opened the modern way of calculating the stopping power of 
swift ions in condensed matter. The case of low- velocity projectiles is much more complicated 
due to the strong interaction of the moving ion with the solid. In this case, the projectile 
is dressed by a number of electrons that strongly screen the ion-solid interaction. Brandt 
H introduced a Thomas-Fermi statistical model in order to define an effective ion charge 
that takes into account how the bound electrons dress the projectile. Other researchers 
|jj have developed the Lindhard approach and have calculated the stopping power using a 
linear-response function. In recent developments |||| , the stopping power for ions at low 
and intermediate velocities has been obtained by introducing the different electron-loss and 
-capture processes associated with the interaction of the projectile with the target [JlOj . 



An important development in the calculation of the stopping power for very-low- velocity 
ions in solids appeared with the application of the local density theory to this field. 
Echenique, Nieminen, and Ritchie (11]] calculated the stopping power for very slow ions 
moving in a uniform electron gas, using well-known techniques in this field. This approach 
has yielded a substantial improvement in the agreement between experimental data and 
theoretical calculations. The main limitation of this approach, as it has been used in the 
actual calculations, is the assumption of having an uniform electron gas in the solid. Al- 



though some attempts are currently tried to improve over this simplification (T^,[T3[, it could 
be convenient to try, at the same time, other alternatives that might be appropriate in the 
case of having very ionic or covalent solids. 

The aim of the work presented in this paper is to apply to the stopping power field an 



approach recently developed for the calculation of electronic properties of solids [14,15 



This is a linear-combination-of-atomic-orbitals (LCAO) approach, whereby the electronic 
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properties of the solids are calculated from the localized wave functions of the atoms of the 
solid. This approach tries to emphasize the local chemical properties of the solid and is 
deeply related to the work done by other groups, trying to calculate the stopping power for 



ions in solids using the stopping power in the vapour target JT6| . The advantage of these 
approaches is related to the non-uniformity of the target, since a local-density-approximation 
(LDA) calculation usually assumes a uniform electron gas inside a crystal. Thus the long- 
term aim of our approach is first, to calculate the stopping power for ions as a function of 
the ion position, in particular, near crystal surfaces; and secondly, to take into account the 
contribution of the different atomic orbitals of the target, mainly those orbitals which have 
such a localized size that can not be replaced by an uniform electron gas. 

In this paper we have chosen to analyze the case of helium interacting with alkali metals. 
This is a case in which the interaction of the projectile and the target is simple. It is, 
however, a complicated system since it presents a long-range interaction between orbitals 
located at large separations. In these metals, the local density approach can be expected to 
be very good; therefore, we have chosen it as a strong test to the method we have developed, 
and the results obtained give strong support to it. 

In the Sees. II and III, we present our model, the formalism used to solve it, and its 
application to the case of He in metals. In Sec. IV we discuss our results, and in Sec. V we 
present our conclusions. 



II. MODEL AND FORMALISM 

A. General formalism 

Our model is an extension of a previous approach to the calculation of the electronic 
properties of the solids using a LCAO method |T^JT5|| . The basic idea is to introduce the 
atomic orbitals ipu, u = i,a, i referring to the crystal site and a to a particular orbital, and 
the orthonormal basis M 
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<t>» = £Os , - 1/2 WVv, (i) 

V 

with 

^ = (^|^), (2) 



obtained using Lowdin's orthonormalization procedure [17]. Using this new basis, the elec 



tron Hamiltonian of a given system can be written in the following way: 

with the operators cl a and c vu related to the ortho normalized wave functions (j) v . The 



different terms in Eq. (|3|) are analyzed in Refs. |14j] and [13] ; here we only comment how 
to introduce the many-body terms of hamiltonian ([3]) in a one-body Hamiltonian by means 
of a Slater-like potential. This implies replacing Eq. (|3]) by the effective Hamiltonian: 

H c z = Y.Kh va + T^(c+ a c va + c+ a c, a ) (4) 

where 

+ E J%(n,-*)+ E (4?-^ t + 4°^i + ^ c )(V) (5) 



V*' c is the exchange and correlation potential |I3| associated with the many-body terms of 
Eq. @. 

We start from Hamiltonian (|j) and assume that its solution can be obtained in the static 
limit v — > 0, for the case of an atom moving inside a crystal (see Fig. 1). In our model, the 
different parameters of Hamiltonian (|]), as well as its static solution, are calculated for a 
geometrical configuration, at each position of the external atom inside the crystal. 

To proceed further, we assume that, due to the atomic motion, there is a time depen- 
dence of hamiltonian (Hf) through the ion velocity. This implies introducing a quasiadiabatic 



Hamiltonian, H e g(i), with the different parameters, E° and T ut± , having an explicit, but 
slowly, time dependence. 

In order to calculate the stopping power at a given time and atomic position, the static 
solution of hamiltonian H g r is introduced. This implies writing 



H ef * | n) = E n | n). 



(6) 



Then, the stopping power (written as a function of the local time t, defining the projectile 

dH c s (t 



position) is given by the following equation |TJ 

-uu„o(t— t ) 



dE 



-2Re£ 



dt'- 



■(o 



dH eS (t) 



dt 



n){n 



dt' 



0). 



(7) 



[We are using atomic units (h = m = e~ = 1) .} Equation (0) is only valid in the quasiadia- 
batic limit, with the ion velocity going to zero. Notice that in Eq. (0), the eigenstates | n) 
correspond to the full Hamiltonian -ff e fr, including the external ion, at the final time t. This 
approximation is obviously only appropriate for v — > 0. 

Equation (0) can be further modified by noting that the dependence of H e g with t 
appears through the coordinate R = Rq + vt, of the external atom. Thus we write dl1 ' " 1/1 



dt 



[v ■ V)-ff e ff(R), and introduce the Fourier-transform iT e ff(q) of if e ff(R). This yields 



dE 
~dt 



<iq da[ 

(2tt) 3 (2tt) 3 

xe iq-(Ro+vt) e -iq'-(R +vt') ( 



-2ReE 



-iw n0 (t-t') 



W n0 



(q- v)(q • v) 



(8) 



] (0\H eff ( q )\n)(n\H eff (c{)\0). 

This equation can be easily integrated on t' . Moreover, we introduce the one-electron eigen 
functions and eigenvalues, | k), of Hamiltonian H e R in Eq. (||) to define | n) and w n Q. 
These steps yield the following results: 

dq f dq' (q • v)(q ' • v) 



— 4 V 

dt ~ V (2tt) 3 7 (2tt) 3 



Wkk' 



k(kp ,k')kp 

x(k' | F cfr (q)e^ R | k)(k | ^ cff (q)e-^- R | k') 
x8{w kk > + q- v), 



(9) 



where the spin has been added up and Wkk' = — £k ■ Note that | k) and are the 
eigenf unctions and eigenvalues of the total Hamiltonian, H e ^{R) , with the external ion 



included. One should remember, however, that | k) and | k') are not eigenfunctions of 
^eff(q) : 

if eff (q) = J dR! e"* R 'H cS (K ')• (10) 

It is of interest to make contact between Eq. ([]) and the linear-response theory. In this 
case, the total Hamiltonian is written as the sum of the unperturbed hamiltonian H$ and a 
perturbation H pcTt = V. Then, Eq. (H) can be transformed by taking for | k) and | k') , the 
eigenfunctions of H$\ moreover, the perturbation V, can be written as follows: 

Z 



V 



where Z is the external ion charge and R its position. Then, the power loss is given by the 
following equation (linear theory): 

dE r dq f dq' (AitZ\ ( AttZ^ 



dt Hk ^ )kF J W J ( 2 -) 3 V ? ) V <z' 2 



iq- v J 



X 



e *(q-q')-R( k ' | p + (q) | k)(k | p(q') | k')S(w kk> + q ■ v) (12) 



dE „ , dq f^\(^\ 



or, equivalent ly, 



dt J (2tt) 3 J \ q 2 J \ q' 2 

x e ^')- R Im X (q,q';-q.v), (13) 

where Inr\;(q, q'; w) is the metal polarizability. 

For an homogeneous system, only q = q' contributes, and Eq. ([i3|) yields 

dE f dq (AnZ\ 2 . XT . 

in agreement with other Refs. [§,0 • Equation (0) is the basic equation giving the stopping 
power of the moving ion, in the low velocity limit, within our LCAO approach. In Eq. (|9|) 
the critical quantity to calculate, using the static interaction between the external charge 
and the solid, is (k | H eS (q) | k'). In this paper we shall concentrate on the He case; this 
provides a simple case in which to test the method discussed here. 
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B. Static interaction of He with a metal 



In this section, we will present a summary of the main results discussed in Ref. ||14|| . We 
shall also extend this discussion in order to calculate the matrix elements (k | if c g(q) | k'), 



needed for the calculation of the stopping power. Following Ref. fTM , we start by considering 



the one-electron interactions between the He Is level and a metal band that is represented 



in Fig. 2 by a half-occupied s level. As discussed in Ref. flljj, there are two different one- 
electron interactions. First, due to the overlap S between the He Is wave function and the 
metal orbital (5* = ("0m I V'He)); there is an increase in the kinetic energy of the electrons of 
the system. This is measured by the following shift of the one-electron terms: 

6E® = 1 -S 2 (E M -E° He )-ST, (15) 

SE^^-^S 2 (E° M -E° ne )-ST, (16) 

where T, the hopping between the two orbitals, ipM and V'He, is found to be — ^S(Ej^ — E# e ). 
E° M and E# e are the metal and He energy levels. Second, due to the hopping T between the 
two orbitals we find a hybridization contribution to the total energy given by the following 
shift in Em and E^ c : 

SE® = (E T \ y (17) 

SEg = - T \ . (18) 
Combining Eqs. ([15| ) -([T8| ) , we find the following contributions: 

5E M = S 2 (E M -E ne ), (19) 

SE Hc = 0. (20) 

These shifts in the one-electron levels yield the following contribution to the repulsive 
energy: 



^K°puisiv e ctr ° n = n M S 2 (E M - E He ), (21) 



where um is the number of electrons in the metal orbital. 



Many-body contributions have also been discussed in Ref. []14 |. These terms can be 
written in a way similar to Eq. (^Tj); in Ref. [14]] it was found that the total repulsive energy 
between the metal atom and He is given by 

^Kepulsi!c Ctr ° n = n MS 2 (E M — -Erc) + n Al{ — Jt + S 2 Jq) + Electrostatic, (22) 

where J° is the exchange integral between the metal and the He orbitals, Jo the Coulomb 
interaction between the same orbitals, and Ki ec trostaric the electrostatic interaction between 
the total charges of the two atoms. For a He-orbital going like (^-) 1,/2 e _/3r we find that 

3 

— J, x + Electrostatic = ~oPS ■ (23) 

O 

This shows that the repulsive potential can be written as follows: 

3 

Eepuisive = n M S (E M — E He — -(3 + J ). (24) 

o 

In our actual problem we are interested in calculating (k 7 | JJ er r(R) | k) , the matrix 
element of the total Hamiltonian between the one-electron states | k) . We will show how 
Eq. (|24"D can be related to (k' | H c s(Tl) \ k) . To this end, we start by discussing the 
solution of the total Hamiltonian (crystal plus the external atom) within a one-electron 
approximation. The solution of this Hamiltonian H is given by 

V> = C k^k + C He 1p R e, (25) 
k 

where ipk are the eigenfunctions of the crystal Hamiltonian, Hq, and ^He the Is orbital 
of He. In writing Eq. (p5|) , we assume that the total Hamiltonian (in our one-electron 
approximation) is given by H = H + Vh 6 > where Vh 6 defines the one-electron potential 
created by the atom. The eigenvalues and the eigenfunctions of H are given by the secular 
equation 



S 



det | (ipi | -E + H | fy) |= 0. (26) 
Now, we follow Ref. [H and introduce the orthonormalized wave functions [as done in 



Eq. (HD for the basis if) u ] 

0i = E( 5r_1/2 )«'^'. ( 27 ) 

i' 

with 

5 fcHe = (V>k I ^He), (28) 

and 

S kk > = (ipk I fa) = 0. (29) 



Using Eq. fl27|) we define the following effective hamiltonian 

# cfr = S- l/2 HS- l/2 . (30) 

In Ref. , the diagonal terms of the effective Hamiltonian were calculated up to second 
order in the overlap, a small parameter used for calculating S 1-1 / 2 in a series expansion, 
while the off-diagonal terms were only obtained up to first order. In our actual problem we 
need to calculate (H e s)kk' up to second order in the overlap, the smallest surviving term of 
the expansion. 

Proceeding in this way, we obtain the following results: 

(-^cff)fcHe = TfcHe = — T^ 1 ^^ ~ ^He); (31) 
= T kk , = {Vu c ) kk > ~ — (Tkn e S}lek' + Wfe*) 

1 E° + E°, 

+ t( ■ En e )S}i e k>Sji e k, (32) 

where E^ and E^ e are the k-state and the atomic levels, respectively. 

Equation was already discussed in Ref. , and found to be valid for a very localized 
wavefunction like the He Is level. Equation (^) is the new equation we are looking for; here 
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(Vffe)fcfc' is associated with the direct perturbation introduced by the He atom on the metal. 
This perturbation is basically due to the atomic Hartree potential, and to the exchange 
perturbation created by the He Is level. 

Equations ([Jl|) and ([32]) can be further approximated by taking E k , the one-electron k 
state levels, equal to E^ a mean level of the metal band (notice that the He level E^ e is 
very deep and that replacing E® by E° M is a good approximation) . Then Eqs. ([H]) and ( 32]) 
read 

TkHe = —-SkHe(E M ~ ^He)> (33) 
Tkk' = (Vllc)kk' — -{TkUeSllek' + ^fc'HeSHek) + J {Em — -^He) ^Hefe'^HeA; 

3 

= (VRe)kk' + T (T/m — E}i c )Siiek'SHek- (34) 



The terms appearing in Eq. fl34|), that depend on T fcHe and SWe, are equivalent to the ones 
going like (—ST) , in Eq. (|I5|) , if T is replaced here by —^S(E M — i?He) j this shows how 
the one-electron correction to the metal level jS 2 (Em — -Ehe) coincides with the one-electron 
contribution to the off-diagonal term in Tkk' if S 2 is replaced by SkHeSuek'- 

Returning to Eq. fl33|), we should comment that T^He is a first order term in the overlap 
Snek while T kk i is second order [(T^He)fefe' included]. The first order term T kiie introduces an 
effective second order contribution to T k y given by 

TkUeTuek' 



(35) 

Combining Eqs. ( p3|) and (|34| ) with Eq. (|35| ) we get the following effective interaction: 

Tkk' = (^HeW + (E° M — E^SkReSnek'- (36) 

This is the one-electron contribution to the effective hopping between the crystal wave 
functions | k) and | k'), as induced by the external atom. When the crystal wavefunctions 
| k) are developed in a local basis 

|k>=5>(k)0 4 , (37) 

i 
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<f>i being the orthonormalized wave functions associated with the metal atom , Eq. fl3~6|) 
reads as follows: 

Tiv = (VHe)ii' + {E M — Eft e )Sm e Sn e i<. (38) 

Equation ( |3"B] ) is the fundamental equation making contact between the repulsive poten- 
tial given by Eq. (|l^) and Ta>. Many-body contributions are partially taken into account in 
Eq. (^) by means of the term (VH e )u' which includes the bare Hartree and bare exchange 
contributions, equivalent to Kiectmstatic and — J° in Eq. (|^). The extra term S 2 Jo, appear- 
ing in Eq. (p2|) is due to the effect of the overlap between the | k) and He orbitals in the 
total exchange interaction. 

This discussion and the results of Eq. (^) suggest to introduce the following effective 
interaction between the i and i' orbitals 

3 

(^eff)^' = SiU e Suei'(E M ~ E Ue — -(3 + (J )). (39) 

This equation should be compared with Eq. fl2l|) that yields the total repulsive potential 
between He and the metal atoms. 

In this equation, (J°) is associated with the effect of the overlap between the He Is 
orbital and the atomic wave functions of the metal in the exchange interaction created by 
the He-orbital. In Eq. ( f24|) , J° is the Coulomb interaction between the atomic wavefunction 
and the He Is orbital; in Eq. ( |3~9"D we have introduced ( J°) , the mean value of this Coulomb 
interaction in the crystal unit cell (the change of J° along this unit cell is small, less than 
10%). 

Equation (|39|) is the main equation giving the effective matrix elements creating the 
excitation between the i and i' orbitals, or the Bloch wave functions | k) and | k') in the 
crystal, in this basis: 

3 

(^eff)ifcfc' = SfcHeShefc' (E M — E Rc — -(5 + { J )). (40) 
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III. DYNAMIC INTERACTION OF HELIUM WITH A METAL 



Once we have obtained the static interaction of He with the metal, and the effective 
matrix elements, we will discuss how to combine this result with the general Eq. (|^) to 
calculate the stopping power for He. First of all, let us mention that we shall use Eq. (|9]) by 
assuming that | k) and | k') are well described, for the He case, by the unperturbed crystal 
wave functions. This is a good approximation in our current case due to the small overlap 
between the He Is and the localized metal wave functions. 

Then, the starting point is the equation 

PeS (R-)]fcfc' = VoSfcHeSfiefc', (41) 

where 

V = (E M -E Rc -l(3+(J )). (42) 

o 

The overlap between the Is He state and the | k) wave functions is written in the following 

way 

(k | ^He) = Jdr ^*(r)^He(r) ~ V£( R He) I dr ^Hc(r), (43) 

where we replace iplir) by ^(R-He) , assuming the He Is level to be very localized. This 
allows us to write: 

(k | H cS (R) | k') = V ^(RHe)Vv(RHe)[/ dv ^He(r)] 2 

= ^(RHc)MR-He). (44) 

This yields [see Eqs. (g), (jig)] 

#cff(q) = VI J rfRe^ R ^(R)^(R) 

= K4*/(q) (45) 

and 
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— = a V (V) 2 f dq - f dq ' v ^ q/ ' v ^ 

dt k{k ^, )kF [ 0) J (2n)i J (2«y Wkk , 



xI kk/ (q)I k M)e^-^ R 5(w kkf + q • v). (46) 

This is the general equation giving the power loss at a given point R. Notice that due to 
the crystal simmetry k' = k — q and q' = q — G, G being a crystal reciprocal vector. 

If we are only interested in the mean power loss and neglect the R dependence, we should 
concentrate on the q = q' contribution. Then Eq. yields 



l-oo (2tt) 3 J-oo (2tt) 3 
X I kk , (q) I k , k (q) S (w kk > + q • v) , (47) 



with k' = k — q and is the step function. 

Equation (|47| ) depends on the velocity direction of the projectile. As we shall only 

consider the case of He moving in alkali metals, crystals that have a very small anisotropy, 

we shall calculate the stopping power by taking an average on all the v directions , which 

will enable us to compare our results with other works fll |. Then: 

1 dE j\dcos9 v dE / dt 
v dx 2 v z 

which redefines dE /dx. 

Equation ([47]) is our fundamental equation for calculating the stopping power for He, 
in the low- velocity limit. This equation can be written in a local basis by developing the 
k states in the atomic orbitals of the crystal. In general, we shall assume that the metal 
wavefunctions are given by an effective one-electron Hamiltonian f/o, such that 

F |k) = £(k)|k>. (49) 

Then, the solution of this hamiltonian yields 

|k) = E c «( fc ) eik ' R, ^( r - R *)> (5°) 

where 4>i Q (r — Rj) are the orthonormalized wave functions associated with the i site (a 
measures the number of orbitals per site). On the other hand ^(r — Rj) should be expressed 
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as a function of the localized atomic orbitals ipi a (r — Rj) using Eq. (1). By substituting 
Eqs. ( |50D and ([I]) into Eq. fl4"?|), we find the following result 

l d l = ^vi? £ £ ^ £ ^ <^e ( * F - * )e( v - m 

x E E <(k)cr(k)4(k')c s (k') 

Ri,R 2 a,/3,7,5 

x E(^(k)- 1/2 )^^(q)(^(k)- 1/2 )^ 

a',/3' 

^(5(kr 1/2 )^S(q)(s(kr 1/2 ) ff 

7', 5' 

^(^fcfc' + q-v), (51) 



xe 



t(k-q)-(Ri-R a 



where 



J#(q) = Jdv e^(r)^ 7 (r - R x ), (52) 
(S(ky 1/2 ) af3 = ^e*' R (S" 1/2 (R)) Q/ 3, (53) 



R 



5(k)^ = 2 e lk ' R / dv ^«(r)^(r - R). (54) 

R J 

Finally, we relate c^(k)c a (k) to the metal Green functions Gp a (\i,w) by the equations 

1 rE F 

®(k F - k)c (k)c* a (k) = - dw ImG^k, w) (55) 

7T J-oo 

1 r°° 

Q{k' - k F )c* s {k)cJk) = — / dw lmG 5 Jk,w). (56) 

7T JEw 



This yields 



1 dE 2 r 1 of 00 dk f°° dq (q ■ v) 



2 



vdx v u/ 7-1 7_oo (2vr) 3 y-oo (2vr) 3 v 

x E E Im[G Q7 (k)]Im[G^ (k)] 

Ri,R 2 a,/3,y,6 

x E(^( k )" 1/2 )-'^(q)(5(k)- 1/2 )^ 

x E(5(kr 1/2 ) 7 y/?Mq)(S(kr 1/2 )«' 

Xe i(k- q ).(R 1 -H 3 ) ( j^ fcfc< + q . v ) ; (57) 
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where 

G/fa(k)= f EF —G Pa (k,w) (58) 

and 

G* r (k')=/ — G^(k>'). (59) 

J oo 7T 

Equation ( |57|) allows us to calculate the stopping power for He in metals, as a function 
of the Green-function components G a p{k) of the metal (calculated in the orthonormalized 

A —1/2 

basis), using a one-electron hamiltonian Hq{k) and the overlap matrix S a/3 (k) associated 
with the atomic wave functions i[) a and ipp . Moreover, also depends on /^(q), the 
Fourier-transform of the overlap between the atomic orbitals ipp^) and ip 1 {r — R) as given 
by Eq. © . 

On the other hand, in order to analyze the stopping power as a function of R we take in 
Eq. (pE6|) q' = q — G, and only the G vectors perpendicular to the v direction. This yields 



for the G component of -4^, 

1 v ax ' 

x4 fc '(q)4a(q - G)e-* G - R 5(wte + q • v), (60) 

and remember that k' = k — q. 

This equation can be written in a way similar to Eq. (f)7]), as a function of S(k) , G a p{k) 
and , 7?3 . For the sake of brevity, we shall only mention here that in general the stopping 
power S = -4^ can be written as follows: 

1 v ax 

5 R = 5 + E^Ge iG ' R , (61) 

G 

as a function of R, where Sq is the mean stopping power given by Eq. (p7[), and Sg the G 
component of Eq. (|60|). Once we have chosen the G vectors perpendicular to v, we have 
calculated Sg by taking an average on the angle between v and q as in Eq. 
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IV. RESULTS AND DISCUSSIONS 



We have applied the previous formalism to the calculation of the stopping power for He 
in alkali metals. For simplicity, the band is assumed to be well described by means of a 



single s orbital. Then, Eq. (51) can be further simplified into the following equation: 



or 



\t - ^ £ /_: ^ /_: ^ - *w - m 

x £ (S(k)- 1 / 2 )/ B -.( q )(S(k)- 1 /')(S(k')- 1/2 )/ R '(q)(S(kr 1/! ) 
Rl,R 2 

xe f(k-q)(Ri-R^( TOM , +q . v ), (62) 



i> ctx J-i J-00 (2n) 6 J-00 {2ti) a v 2 

Rl,R.2 \ R 1 / \ R 2 

x5(w fcfe / + q- v), (63) 



where 



and 



J R i(q) = f dr e iq r ^(r)^(r - Ri) (64) 



S(k) = J2e ikR jdr V(r)V»(r-R). (65) 



R 



Equation (|64j) can be written in a more symmetric way as follows. Take 



J R i(q) = e ^- R i/2y rfr e iq -( r - Rl / 2 V(r)^(r - 



3 iq-Ri/2 



/ Rl (q), (66) 



then 



£ e i(k-q).R 17 R l(q) = ^ e i(k- q /2). Rl j R l(q) (6?) 
Ri Ri 
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and 

^ e - l( k- q) .R 2/ R 2(q) = ^ e - l( k- q /2).R 2 J R 2(q)> (6g) 
R 2 R 2 

Equations (|62[) and ( |63j) yield the stopping power for He as a function of S(k) and 
J R (q). S(k) has been calculated using the atomic wave functions given in Ref. ]19j . 
The calculation of J R (q) is more complicated since the Fouriertransform of the atomic 
wavefuncions centered on different sites are needed. This is the well-known problem of mul- 



ticenter integrals. Several solutions have been tried in the literature |2(| such as expanding 
the Slater-type basis in a Gaussian one ]2~T| , |2"2| . We have used, however, an adaptative 
algorithm of integration by Monte Carlo techniques |23j to perform / R . 



An approximate solution, which yields good results for So, is obtained by replacing 

/ R (q) ~ 5(R) * /(q), (69) 



where 



and 



S(R) = Jdr ^(r)^(r - R) (70) 



/(q) = Jdr e iq - r ^(r)^(r). (71) 



Equation ( |S9"D is exact in the limit Rj = or q — > 0. In general, we expect Eq. QBTS| ) to 
give a good approximation to J Ri (q) if q • Rj/2 is small. 
Introducing Eq. d69|) into Eq. (|63D yields: 

x - d 4- = 2,(y>y [ l dcos9 v r*r«* ^le(k F - k)e(k> - k F ) 



vdx v u W-i "J-oo (2tt) 3 7-oo (2tt) 3 V 2 

x W)M) 1 ^ 5(R)ei(k ~ q/2) ' R |2 ^ + q • v) - (?2) 

Equation ( |72|) is the basis of our approximation to Eqs. fl6"2|) and (p3|). We should also 
mention that Ek (the electron energy band of the alkali metal) has been assumed to follow 
a free electron dispersion law. 
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Before discussing the numerical results given by Eq. (|72|), it is worth considering the 
results obtained by neglecting all the overlaps between the alkali atom wavef unctions. Then 
we write 



S{R) = { ' (73) 




and 

S(k) = 1 (74) 
and replace Eq. ( |7^ ) by the following equation: 

_ = 2tt (V 2 / dcos6 v / — - / — \ (k F - k)G(k' - fc F 



X 



J dv ^ 2 (r)e iq - r2 5(^ + q-v). (75) 



It is also convenient to discuss at this point the stopping power given by the following 
simple model: a uniform electron gas interacting with a slowly moving He atom by means 
of the following contact potential 

#pert = W-Vt). (76) 



Here V is assumed to be the same local potential introduced in Eq. (f44|). It is an easy task 
to develop this model following the same steps as discussed above for the LCAO approach 
and find the following expression for the stopping power 

\t = ^ L w? jC wf { ^ e{kF - k)e(k ' - k ^ + q • v) - (77) 

The integral in cos9 v equals 2 because in the latter expression depends only on | v | . 

Comparing Eqs. ( [75D and (|77|) , we see that their only difference is associated with the 
term | I(q) | 2 =| / dr ip 2 (r)e l(l ' r | 2 , which gives the form factor of the metal orbital. We 
should also comment, regarding Eq. (|72|) that in the alkali metals S(k) ~ S(kp) since, in 
the low velocity limit we are considering k and k' are located near the Fermi sphere, S(k) 



being almost constant on this surface that presents a very small anisotropy. Then, Eq. (|72|) 
can be obtained from Eq. (|75| ) by replacing the form factor I(q) by 

where 

S(k - q/2) = J2 5(R)e i(k ^ q/2) R . (79) 

R 

Thus the three different cases we are considering yield the same equation for the stopping 
power, but for a specific factor taking the values 1, 1(q) and -D(q, k), for the free-electron gas 
(FEG), the LCAO model with S(R) = for R ^ (LCAO-I), and the LCAO model with 
S'(R) 7^ (LCAO-II), respectively What is of interest to realize about this discussion is 
that the free-electron-gas model overestimates the stopping power, while the simplest LCAO 
model underestimates it. In Table p], we give the three values of the mean stopping power So, 
for He in Na as calculated from these equations. As shown in this Table | the free electron 
gas model yields a stopping power three times too large, while in the LCAO-I model is 
about eight times too small. 

One word of caution must be put here. The FEG model discussed here can not be 
compared directly with the LDA used to calculate the stopping power of He in metals. The 
point to notice is that in the model defined by Eq. ([T6|), V ' is the contact potential for the 
interaction of He with the s orbitals of the alkali-metal atoms. The model of Eq. ([T6|) is 
only introduced here in order to explain how the form factor of Eqs. ([78|) or (|75| ) is the main 
term controlling the He stopping power. 

As regards the factor -D(q, k) used to calculate in the LCAO-II approximation, 
notice the strong dependence that -D(q, k) has on the number of neighbors used to calculate 
S(k) = E R S(R)e ik ' R and S(k - q/2) in Eq. fl78|). We have found that in order to get 
a reasonable accuracy (around 5%) it is necessary to add up to the fifth or sixth neighbor, 
depending on the alkali metal. 



As mentioned above, Eq. ( |SBD has been accurately calculated for Na using Monte Carlo 
techniques. We have found that this Monte Carlo calculation yields 
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1 dE 

— — = 0.085 a.u. (Na), (80) 
v ax 

a value a little larger than the one found using our LCAO-II approximation. By assuming 
the same correction factor for all the alkali metals, we find the results given in Table |J 
column (a). This Table also shows the theoretical figures obtained by Echenique, Nieminen, 
and Ritchie |11 . 



We see from Table O, column (a), that the results for K and Rb are in excellent agree- 



ment with Ref. fil l) although the stopping powers we find for Li and Na are a little larger. 
This difference can be partially attributed to the simple model we are using, since a single 
s orbital per alkali-metal atom has been assumed to form the metal conduction band. This 
approximation can be expected to be a reasonable one for very electropositive atoms like K 
and Rb , but not so appropriate at least for Li. Thus, in the calculations of Papaconstan- 
topoulos [24|] for Li, only 52% of the occupied density of states has a s like character. If we 



introduce in the results of Table II, a factor 



<l< (Rb) (81) 

which normalizes the stopping power of each alkali metal to the total number of s electrons 
with respect to Rb, we find the results of Table 0, column (b), in much better agreement 
with the LDA calculations. 

The conclusion we can draw from these results is that the method developed in this 
paper is quite appropriate to calculate the stopping power for He moving slowly in alkali 
metals. We can also expect that the method will be useful to calculate stopping powers for 
atoms in transition metals. 

In a further step we have calculated, using Monte Carlo techniques, the stopping power 
dependence on the ion position (for He moving in a channeled direction). We have considered 
that He moves in a Na crystal along the [100] direction. We have calculated the different G 
reciprocal vectors contributing to the stopping power [Eq. (|60D 1; this implies taking the G 
vectors perpendicular to the [100] direction. In a bcc lattice, the first reciprocal vectors to 
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be considered are the followings: G = ^(0, 1, 1), ^(0,0,2) etc. Using Eqs. (0) and (P) 
we have obtained the stopping power Fourier components So shown in Table [HI]. 

Figure 3 shows Sr., with R changing in a [100] plane. The main conclusion we can draw 
from these calculations is the strong dependence that the stopping power shows as a function 
of the impact parameter: the stopping power can vary as much as 100% for different impact 
distances. We should comment that these changes are not associated with the electronic 



metal charge [|13l ; this charge, as obtained in our LCAO approach with an s level per atom, 
appears to be almost constant in the crystal lattice except very close to the atomic sites. 
For a He-atom channeled along the Na [100] direction, one expects some kind of oscillatory 
motion of the atom, with the impact parameter changing along the He trajectory. Then, 
the mean-stopping power for the channeled case would appear as an average of the different 
values shown in Fig. 3 around the minimum value of the stopping power. Each case should 
be analyzed specifically, but assuming the incoming atom to explore only half of the total 
available space, one would get around 50-60% of So, namely 0.04 a.u., 80% of the value 
calculated in LDA. 

V. CONCLUSIONS 

The aim of this work has been to develop a first-principles, free-parameter, approach 
based on a LCAO method to calculate the stopping power for atoms moving in condensed 
matter. In the past few years the interest in, generally speaking, tight-binding methods 
|25| for analyzing the electronic properties of solids has increased a lot. This emphasis is 



partially due to the interest in using a local point of view, closely related to the chemistry of 
the local environment. The work presented in this paper follows this general trend and tries 



to apply the ideas recently developed in Refs. p^ , [15|j for analyzing the electronic properties 



of solids following a LCAO method, to the stopping power area. In the long term, this 
approach can be expected to be also useful for analyzing other dynamical processes like the 
charge transfer between moving ions and the solid, sticking mechanisms, etc. 
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In Sec. II, we have presented our general approach and have related the stopping power 
for atoms, in the low- velocity limit, to the electronic properties of the crystal as described 
using a LCAO method. All the parameters appearing in Eq. (^), the general equation giving 
the stopping power, can be obtained from the local wave funcions of the atoms forming the 
crystal. Equation (|Sp has been applied to the case of He moving in alkali metals. He is 
a simple atom, but the alkali metals present a strong test to our method as their atomic 
wave funcions interact strongly with each other up to large separations. In Sec. IV, we have 
presented our results and have found that the stopping power for He is very well described 
with our local LCAO approach, if the interaction between different alkali metal atomic 
orbitals is included, at least, up to fifth neighbor. 

We conclude that the LCAO method discussed in this paper offers the possibility of 
calculating accurately the stopping power for ions moving in solids. This could be a con- 
venient framework for analyzing solids having localized d bands and for discussing specific 
geometries like the case of atoms moving near surfaces, or the channeled case discussed in 
Sec. V. 
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FIGURES 



FIG. 1. He moving inside a metal. 

FIG. 2. Schematic representation of a He atom interacting with a metal band simulated by a 
metal level Em- 

FIG. 3. Stopping power for He moving inside a Na crystal along the [100]-direction. The 
stopping power is normalized with respect to its mean value. The coordinates correspond to a 
(100) plane perpendicular to the projectile trajectory. The atomic rows along the [100] direction 
are projected onto the points having the coordinates: (0,0), (0,1), (1,1), (1,0) and, (0.5,0.5). 
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TABLES 





FEG 


LCAO-I 


LCAO-II 




0.140 


0.007 


0.056 



TABLE I. Stopping power for He in Na as calculated for the free electron model (FEG), the 
LCAO model neglecting the metal wave functions overlaps (LCAO-I), and the LCAO model taking 
into account these overlaps (LCAO-II). 



(a) 



(b) 



(c) 



Li 
Na 
K 
Rb 



0.260 
0.085 
0.026 
0.015 



0.110 
0.068 
0.023 
0.015 



0.100 
0.053 
0.023 
0.016 



TABLE II. Stopping power in a.u. for He in different alkali metals, (a) Our results, as 
calculated using Monte Carlo techniques, (b) Our results with the s ocupancy correction, (c) 
ENR, from Ref. p|. 



atomic units 



S(0,0,0) 0.85 x 10" 1 

S(0,1,1) 0.18 x 10- 1 

S(0,0,2) -0.04 x 10- 1 

TABLE III. Stopping Power for He in Na for the channeled direction [100] 
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